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Abstract 

' We propose to quantify the "correlation" inherent in a many-electron (or many- 

fermion) wavefunction ip by comparing it to the unique uncorrelated state that has the 
^ ! same 1-particle density operator as does 

m 

^ ! Electron correlation is of fundamental importance in quantum chemistry and magnetism, but 
as yet there is no definitive way to quantify it: given the wavefunction ip(xi, x 2 , ■ ■ ■ , x^) repre- 
senting the state of a system of N electrons, how much "correlation" is there in that iV-electron 
state? By definition, a wavefunction that has the form of a Slater determinant represents an 
"uncorrelated" state, but how much "correlation" should be attributed to states that are not 
represented by a Slater determinant wavefunction? Some measures of "correlation" have al- 
qh ready been advanced in the literature: the "degree of correlation" |TJ |2| and the "correlation 



entropy" j3J EJ El EH El , for example. These correlation measures depend completely upon the 
eigenvalues of the 1-particle statistical operator 7, i.e., the operator with the integral kernel 



7(x, y) = N j ip(x, z 2 ,..., z N )ip(y, z 2 , ■ ■ ■ , z N ) dz 2 dz 3 ■ ■ ■ dz N (1) 

(supposing that the wavefunction ip has norm 1, so that the trace of 7 equals the number of 
electrons). The "degree of correlation" is inversely proportional to the sum of the squares of 
the eigenvalues of 7, and the "correlation entropy" is the Shannon entropy of those eigenvalues. 
Such measures ascribe the same amount of "correlation" to all wavefunctions that have the same 
1-particle statistical operator. We feel it imposes a severe conceptual limitation to have to say 
that all wavefunctions having the same 1-particle operator (or the same 2-particle operator, 
for that matter) contain the same amount of correlation. Here we propose a new measure of 
correlation which does not suffer that limitation. 
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States represented by Slater determinant wavefunctions are the only pure states are deemed 
"uncorrelated" ; but certain mixed states should also be regarded as "uncorrelated" , namely, 
the mixed "quasifree states." These have random particle number and must be represented by 
density operators (i.e., statistical operators of trace 1) on the fermion Fock space. Let Tï be 
the 1-electron Hilbert space, and let a\ and a/ denote the electron creation and annihilation 
operators for / 6 Ji. A density operator F on the fermion Fock space over 7i represents a 
"quasifree state" jS] if 

Tr ( r 4i a /2 ' ' ' a U a 9n ■ ■ ■ a 92 a gí ) = 6 mn det (Tr(Ta\a g] )) n i j=1 (2) 

for all n and all fi, ■ ■ ■ , f m , gi, ■ ■ ■ , g n £ W (6 mn denotes Kronecker's delta). The two-point 
correlations TrÇTa^ag) determine all higher correlations under the quasifree state T, and in this 
sense a density operator satisfying @ may be called "uncorrelated." The 1-particle operator 7 
defined by 

(0,7/) = Tr(Ta\a g ) (3) 

has trace equal to the average particle number under T, i.e., Tr(7) = Tr(rN) where N denotes 
the number operator. 

In case the average particle number is finite, we can reconstruct T from the 1-particle operator 
7 as follows. Let {(fii}^ be a complete system of eigenvectors of 7 with 

7 (0i) = Xí4>í . (4) 

The eigenvalues Aj are all between and 1; we will intèrpret them as "occupation probabilities." 
Let SS denote the class of all finite sets of positive integers, including the empty set. We can 
choose a random member of SS, that is, we can form a random set s of positive integers, by 
including i in s with probability Àj (and excluding it with probability 1 — Àj) independently of 
all other positive integers which may or may not be included in s. This procedure produces the 
set s with probability 

p(s) = n^nc 1 -*) • (5) 

«es -i^s 

For any nonempty s = {ii, . . . , i n } G SS, let \l/ s denote a normalized Slater determinant in the 
natural orbitals 4>n , ■ ■ ■ > 4>i n (the phase of \l/ s is arbitrary) and define 

n— particle space 
Ps = 0© •■■©()© |^s)(^s| ©0©0©··· 

- a projector on the Fock space; in case s is the empty set, let P s denote projection onto the 
vacuum space. It can be shown [ü] that 

r = J>(s)p s (6) 

s£SS 
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satisfies (J2J) and (jïïj). Slater determinants are special cases of quasifree states: if \1> is a Slater 
determinant in the 1-electron orbitals ipi, . . . , ipw then the projector 

N—particle space 

p s = o©···©o© |*)(*| ©oeo©--- 

represents a quasifree state whose 1-particle statistical operator is the projector onto the span 
of {V>i,...,tM- 

Our idea is to quantify the "correlation" in a many-electron wavefunction tjj by comparing 
it to the unique uncorrelated (quasifree) state V having the same 1-particle statistical operator 
7 as does The more \ip)(ip\ resembles T, the less "correlated" we consider tp to be: 

we would like to quantify the correlation inherent in tjj by some measure of the dissimilarity 
of and T. There are several ways one might measure this "dissimilarity" and thereby 

quantify electron correlation; our choice is the following one: we identify ip with an iV-particle 
vector in the Fock space and define 

Corr(^) = -log (ip,Tip) , (7) 

where T is the quasifree density (JUJ) determined by the 1-particle statistical operator (JTJ) for ijj. 
This quantity is nonnegative, and it equals if and only if ip is a Slater determinant. It can be 
seen from the discussion in the following paragraph that Corr(^) < oc. 

It is easy to calculate Corrf^) if the expansion of ip in Slater determinants in the eigenvectors 
of the 1-particle statistical operator (0) is available. The eigenvectors of 7 are called natural 
orbitals and the corresponding eigenvalues lie between and 1 (TO]. Here, "natural orbitals" 
will refer only to eigenvectors of 7 with nonzero eigenvalues. Recalling the notation in formulas 
(HJ) and (JüJ) above, let SS^ consist of all sets of N positive integers, each of which is the 
index of a natural orbital. The wavefunction ip is a superposition of Slater determinants in 
the natural orbitals, even when the natural orbitals don't span the whole 1-particle space ^T] . 
Thus ip = J2sçss N c ( s )^s , where \I/ S denotes a Slater determinant in the orbitals indexed by 
s, and we have 

Corr(V>) = -log P( S M S )\ 2 > ( 8 ) 

sess N 

where p(s) is as defined in (jSJ). 

Remark 1: a formula for Corr(í/>) in 2-particle systems. 

When ip(x,y) = —ip(y,x) is a 2-electron wavefunction Coty(i/j) is a functional of the eigen- 
values of the 1-particle statistical operator. In this case, it is known ^2] that there exist 
an ortho normal system <7i, /2 3 <72> • • • } of 1-electron orbitals and nonnegative numbers 
pi,P2, . . . such that 

ip{x,y) = VPj^(fj( x )9j(y) - 9j(x)fj(y)) ■ (9) 

3=1,2,- 
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The wavefunctions fi and are natural orbitals with occupation probabilities pi, for 

7 = Y^Pi(\fi){fi\ + \9i){ffi\) ■ 

i 

In this case, one can calculate that 

Corr(^) = -log Y^Pi{pil[0--Pj)Y ( 10 ) 

using formula (JSjl. Formula (jl(J|) shows that Corr(^) can be arbitrarily large, for the 2-electron 
wavefunction (jïïj) can have arbitrarily small coefficients pj ^3] • 

It is interesting to apply formula (|l(Jj) to an exactly solvable model: the "two-site Hubbard 
model" that was investigated in [ü] as a study of "correlation entropy." This simple model 
can be used to illustrate the "kinetic exchange" correction to the Heitler-London theory of 
diatomic hydrogen [T3j. There are two fixed spatial orbitals, labelled 1 and 2, whose span must 
accomodate two electrons, and the Hamiltonian is 

H = —t (a\^a 2 ï + c4t a iT + a \i a n + a 2| a u) + U ^ a ít a iT a í| a U + a 2T°2T a 2| a 2J.) > 

where U is the on-site repulsion energy (or attraction energy if U < 0). The 2-electron ground 
state of H depends only upon the dimensionless interaction parameter u = U/t. Denote this 
ground state by ip u . The ground state ipo for u = is a Slater determinant, and ip u tends 
towards the maximally correlated Heitler-London state ■^( a \-\ a 2i + a 2t a í|)l ) as u tends to 
+oo. The correlation entropy (i.e., the Shannon entropy of the spectrum of the 1-particle 
operator) of ijj u was found [ü] to be strictly increasing in \u\, approaching its maximum possible 
value as \u\ — >■ oo. The same is true of Cott(i/j u ). However, Cott(i/j u ) seems better behaved 
than the correlation entropy S(i/j u ) for small vàlues of \u\: the former is infmitely differentiable 
at u = but the latter is only differentiable once there, and hm S(ip u )/CoTr(ip u ) = oo. If 

S(ip u ) is normalized so that lim S(ip u )/Corr(ip u ) = 1, then S(ip u ) is larger than Corr(ip u ) for 

all \u\ > 0. 

Remark 2: Corr(^) is not a function of the 1-particle operator. 

We have seen that Corr(?/>) is a function of the spectrum of the 1-particle statistical opera- 
tor in case ip is a 2-electron wavefunction. In general, however, Corr(-?/>) is not a function of 
the 1-particle statistical operator. To show this, we exhibit two 3-particle wavefunctions that 
have the same 1-particle statistical operator but contain different amounts of correlation. Let 
ei,e 2 , • • • ,e 6 be six orthonormal 1-electron orbitals, and let us denote 3-particle Slater deter- 
minants in these orbitals by listing the three indices involved (in increasing order) in between 
vertical lines, so that, for example, 

|245| = — =(|e 2 e 4 e 5 ) + |e 5 e 2 e 4 ) + |e 4 e 5 e 2 ) - |e 4 e 2 e 5 ) - |e 2 e 5 e 4 ) - |e 5 e 4 e 2 )) . 
V6 
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Consider the two 3-particle wavefunctions 



| |135| + 



| |246| 



§ (|123| + |345| + |156|) 



These have the same 1-particle statistical operator 
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but Corr(í') w 4.08 and Corr($) « 5.51, taking logarithms to the base 2 in formula (|8*|). (The 
matrix above displays only the matrix elements (e^e.,); outside this finite block all matrix 
elements are 0.) 

Similar examples can be contrived to show that Corr is not a function of the m-particle 
statistical operator for any m > 1 either [T3j . 

Remark 3: Corr for mixed states. 

The measure (J7J) of fermion correlation can be extended to from pure states (given by wave- 
functions) to mixed states (given by density operators). If D is a density operator, we would 
require "the correlation in D" to equal if and only if D is a quasifree density, and of course 
"the correlation in D" must equal Corr(^) in case D = These requirements are met by 



Corr {D) = -2 log (Tr(D 1/2 I\D 1/2 ) 1/2 ) 



where T is the quasifree density (J2J) corresponding to the 1-particle statistical operator 7 defined 
by (0,7/) = Tr(Da\a g ) [TÏÏ]. Formula (HU) uses the generalized "transition probability" ^7] or 
"fidelity" [H] between D and V: 



Tï(D 1/2 TD 1/2 ) 1/2 = \\D 1/2 T 1/2 



I trace 



^1/2/^1/2)1 



/2 



This generalized transition probability recommends itself as a measure of the closeness of any 
two density operators — D and T in this case — because it enjoys several properties that 
distinguish it as a useful quantity in quantum information geometry [Tïïl IT71 ÜU[ l2"Tj. 

Note that (fTTj) assigns positive correlation to most mixtures of uncorrelated states. For 
example, if <pi and 02 are orthogonal 1-particle wavefunctions, then 



Ccyrr(\\4>x){<j>x\+\\(j>2){(j>2\) 



5 



(taking logarithms to the base 2). It may seem odd to assign positive correlation to a 1-particle 
state, but it is due to our identification of "uncorrelated" with "quasifree." This means that a 
state is deemed "uncorrelated" only if the occupations of its natural orbitals are statistically 
independent. In the example, the natural orbitals <f>\ and <p2 are each occupied with probability 
1/2, but the probability that both natural orbitals are occupied (at the same time) is not 
1/2 x 1/2. We consider the state to be "correlated" because the occupation numbers of its 
natural orbitals are correlated random variables. 

Conclusion 

We have proposed a new way to measure the "correlation" in states of many-electron (or 
other many-fermion) systems. The idea is to compare the state to its unique reference state: 
the quasifree state with the same 1-particle statistical operator. To quantify the agreement 
between the many-electron state and its reference state, we have chosen the negative logarithm 
of the overlap, viz. formulas (J7j) and (jllj) . 

Our correlation measure, as formulated above, applies only to finíte systems (àtoms and 
molècules) but not infinite systems (sòlids). Other interesting measures of correlation are 
available for (infinite) uniform electron gases |3J 1221 H3] • Can a serviceable formulation of our 
correlation measure be found for infinite systems as well? 

The most intriguing property of the measure Corr(^) is the fact that it is not a function of 
the 1-particle statistical operator: wavefunctions with the same 1-particle statistical operator 
may have different amounts of "correlation." For wavefunctions of 2-electron systems, however, 
Corr(^) is a function of the eigenvalues of the 1-particle operator, viz. formula (|10|). Further 
investigation of the comportment of Corr(^) is underway, to see whether it conforms to what 
we would intuitively expect from a measure of electron correlation. 
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